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Abstract 

We study the phenomenon of additional light waves (ALWs), observed in crystal 
optics: two or more electromagnetic waves with the same polarization, but different 
refractive index, propagate simultaneously in a isotropic medium. We show that 
ALWs are common in relativistic hydrodynamics, and in particular in strongly cou- 
pled systems that admit a dual gravitational description, where the ALWs are dual 
to quasi normal modes in the AdS gravity. We study both the transverse and the 
longitudinal light wave propagation. In the longitudinal channel we find a transition 
between regimes with different number of excitonic resonances which resembles the 
transition to standard optics observed in crystals. 
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Introduction 

In the recent years an intense research on light propagation in some artificial media, called 
metamaterials, has been developed pQ. Some very uncommon properties, theoretically 
predicted, have become viable. 

One of the most attractive is the negative refraction of light, originally predicted in [2]. 
It is the phenomenon for which the phase velocity and energy flux of a electromagnetic 
wave in a material have opposite directions of propagation. One may ask if there exist any 
generic class of media, for which the EM response functions are exactly computable, that 
have a negative refractive index. With this motivation, in [3] , we studied negative refrac- 
tion in relativistic hydrodynamics and in particular in strongly coupled relativistic media 
We established a new connection between optics and geometry through the gauge/gravity 
correspondence [4j, and we observed that negative refraction is quite ubiquitous for this 
class of media. The negative refraction of light appears because the dispersion relations 
are dissipative and they violate the locality, i.e. there is spatial dispersion. 

Spatial dispersion is an interesting property of media that leads to new amazing non- 
standard optics. In media with spatial dispersion there is a linear non local relation 
between the electric field and the electric induction. This translates in the dependence of 
the permittivity on the wave vector. 
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A non standard optic phenomenon associated to spatial dispersion is the propagation 
of additional light waves (ALWs) in crystal optics. The existence of ALWs has been 
predicted in 1957 by Pekar [5J. He claimed that spatial dispersion effects are relevant 
in presence of resonances given by electron-hole bound states, called excitons. At the 
resonance frequency the interaction among the excitons and the photons is strong, and the 
light wave that propagates in the material is a mixture of these states, called polaritons. 
The energy necessary to create the resonance is proportional to the momentum of the 
incoming photon. As a consequence the linear relation between the electric field and the 
electric induction violates locality, and the permittivity has a pole that depends on the 
wave vector. 

In absence of strong spatial dispersion the dispersion equation has a single solution 
for n 2 , the square of the refractive index. If there are strong spatial dispersion effects 
new solutions to the dispersion relation arise. At the resonance frequency these new 
solutions are associated to the propagation of ALWs. This exotic picture, in which many 
different propagating light waves are associated to the same incoming photon, has been 
subsequently observed in experiments (see [6] and references therein). 

The lifetime of an exciton depends on the order parameters of the theory. If its 
lifetime becomes too short, the ALW associated to this resonance does not propagate 
anymore. Pekar's theory predicts that by varying the order parameters the number of 
waves propagating in the medium varies too. This transition among different optical 
regimes has been observed in crystals with one exciton. In that case for low temperature 
the lifetime of the exciton is large enough for both waves to propagate. As the temperature 
increases the exciton lifetime reduces and above a critical value only the standard single 
wave propagates. 

In this paper we investigate the existence of ALWs in relativistic hydrodynamics and 
in particular in strongly coupled media with a gravity dual description. In such media the 
propagation of light, in the hydrodynamical regime, typically shows spatial dispersion. We 
revisit the model studied in [3], where the refraction of transverse light waves was studied. 
It was found that at finite temperature and chemical potential the refractive index was 
negative in the low frequency and small wave vector-regime. The refractive index was 
computed by considering only the lowest order dependence of the permittivity from the 
wave-vector. Here we observe that if we take into account the whole spatial dispersion 
effects there is a second solution of the dispersion equation, as in the case of crystals. We 
study this additional light wave and we show that it propagates for every value of the 
order parameters. 

Then we analyze the propagation of longitudinal electromagnetic waves in the same 
system. In absence of spatial dispersion longitudinal light waves do not propagate in 
a medium. On the contrary here the permittivity is strongly dependent on the wave- 
vector, and there is a propagating wave at finite temperature and zero charge density. If 
a chemical potential is turned on, there is another solution to the dispersion equation, 
and hence an ALW. The analysis of the propagation of these two waves depends on the 
values of the order parameters. Indeed, by varying them, there is an optical transition, 
between the regimes with one and two propagating waves. This is similar to the transition 
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to standard optics predicted by Pekar's theory. 

The paper is organized as follows. In section [1] we review the linear response theory of 
media with spatial dispersion. In section [2] we discuss the propagation of the light waves 
in media with spatial dispersion and we review the basic properties of Pekar's theory. In 
section [3] we observe that if the medium is described by relativistic hydrodynamics, then 
ALWs are generated. In section H] we study the ALWs in an holographic theory in the 
hydrodynamical regime. In section [5] we study some properties of these waves and then 
we conclude. In the Appendix [A] we discuss the relation between the current correlators 
and the linear response function of a medium with spatial dispersion. In the Appendix [B] 
we discuss the Poynting vector and its relation with the sign of the refractive index for 
the ALWs. 



1 Linear response and spatial dispersion 

In this section we review the linear response electromagnetic functions for media with 
spatial dispersion [7HTU] i.e. we take into account the momentum dependence. In this 
case it is common to use the E, B, D approach to macroscopic electromagnetism for which 
the Maxwell equations in the Fourier space are 

k-B = , k-D = , k A E = ujB , k A B = -uD (1) 

and Di = €ij(u, k)Ej. The linear response of the medium is encoded in the permittiv- 
ity tensor, €ij(u,k). For isotropic media this tensor naturally decomposes in transverse 
€t(u>, k) and longitudinal part €l(u>, k), and 

A = etjiu, k)Ej = e T (oo, k)Ej + e L (u, k)E^ (2) 

where Ef and Ef are the transverse and longitudinal components of the electric field, 
and they are defined in terms of the projectors P? and P-j as 

B? = PSE j =(6«-^E j , Et = P$E j =*g-E j (3) 

The functions er(u),k) and €i,(u>,k) are obtained in linear response theory by applying 
an external electromagnetic field potential A, and measuring the electromagnetic current 
Jf. Ji = q 2 GijAj, where the Gij is the retarded current correlator in the medium and q 
is the four-dimensional EM coupling. As we show in the Appendix [A] the transverse and 
longitudinal permittivity are related to the transverse and longitudinal retarded current 
correlator by 

4:71 47T 

e T (u,k) = l- — q 2 G T (u,k) , e L (u,k) = l- — q 2 G L (u,k) (4) 

Thermodynamical stability requires the positivity of the outgoing heat flow. This is 
associated to the quantities 

Qheat = ^-Im(e T (w, k))\E T \ 2 , Q hcat = ^-lm(e L (u, k))\E L \ 2 (5) 
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for transverse and longitudinal waves respectively. They are positive because the func- 
tion €t and ei are computed from the retarded Green functions as in PJ and both 
Im(GV(w, k)) < and Im(Gi(a;, k)) < 0, which implies Im(eT(u), k)) > and Im(eL(o;, k)) > 
0. 

Using Maxwell equation (0Q), the dispersion relations for the transverse and longitudinal 
waves read 

k 2 

e T (uj,k) = — , e L (u,k) = (6) 

Due to spatial dispersion €t(w, k) and el(uj, k) are functions of both u and fc. Differently 
from [3], here we do not expand at k 2 order. In such a way we can capture the effects of 
large spatial dispersion, i.e. we allow for large momentum dependence in the permittivity. 
By using the dispersion relation ([6]) one can find the on-shell ratio k 2 /u 2 and from this 
ratio the refractive index n 2 (u). This is in general a complex quantity that takes into 
account the dissipation of the medium. Note that, since the permittivity depends on the 
wave vector, the dispersion relation can give in general more than one solution for n 2 (u). 

For a propagating electromagnetic wave, the sign of the refraction is given by the 
reciprocal sign of the phase velocity (Re[n]) and of the flux of energy (the Poynting 
vector). In presence of strong spatial dispersion and strong dissipation, the definition of 
the Poynting vector and the determination of the flux of energy are not well established [7] . 
Hence in this case the sign of the refractive index cannot be study without the knowledge 
of a microscopic description of the system. We discuss about this issue in the Appendix 

m 



2 Pekar's theory and ALW 

In 1957 Pekar [5] discussed the existence of an additional light wave in the optical spectra 
of crystals in presence of strong spatial dispersion. This phenomenon is associated to the 
presence of a resonance in the crystal, namely an exciton. See [UJ for review. For example 
in a semiconductor crystal, an exciton is associated to the transition of an electron, from 
the valence to the conduction band. When such a transition occurs a positive charge 
carrier, a hole, remains in the valence band. This hole and the electron interact throughout 
a Coulomb potential. The bound state generated by this interaction is the exciton. This 
exciton couples to light and its coupling can be described at macroscopic level through 
the dielectric tensor e(u, k). 

The exciton has the same momentum of the absorbed photon, i.e. the energy necessary 
to create the exciton depends on the wave-vector k. This translates in the k dependence 
of the dielectric tensor. In the simplest case, a single exciton, with effective mass M and 
inverse lifetime T, modifies the transverse conductivity as 



e T (w, fc) = e H — — (7) 



A_ 

u o + m ~ u ~ iT 

where uq is the energy required to create an exciton without momentum. 
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On shell the dispersion equation is &r(w, k) = K. Solving this equation and using the 
definition n 2 = K for the refractive index one finds 



w 2 e + 2M(w + iT - co ) ± \ 8AMw 2 + (w 2 e - 2M(w + iT - u )) 2 
•*> = (■» 

In [5] Pekar concluded that the two solutions of the dispersion relation in presence of an 
exciton are associated to two light waves, identically polarized. If there are multi-excitons 
resonances (other poles depending on the wave-vector in the dispersion relation) there is 
an ALW associated to each pole. These waves are distinguished by their different velocity 
and damping. This phenomenon is similar to the birefringence for anisotropic crystals. 
The main difference is that in the case of birefringence the two waves with different 
refractive index have also different polarization. Moreover birefringence is strongly related 
to anisotropy, while ALWs appear in isotropic systems too. 

ALWs can be studied also in longitudinal waves. In that case, in absence of spatial 
dispersion, there is no light propagating in the crystal. Otherwise if spatial dispersion is 
taken into account the propagation is possible, and every pole involving k in ei(u,k) is 
associated to a light wave. 

Propagation of ALW 

The propagation of the ALWs in the crystal does not directly follow from having multiple 
solutions to the dispersion relation. Indeed at frequency far from the excitonic resonance 
the photon does not mix with the exciton, and a single light wave propagates. The 
problem of propagation of the ALW is related to the choice of the boundary conditions. 
Indeed in the case of strong spatial dispersion the usual Maxwell boundary conditions are 
not enough to determine the amplitudes of the reflected and transmitted waves. Pekar 
faced the problem of propagation of the ALWs by introducing an additional boundary 
condition (ABC). There are many possible phenomenological choices of these boundary 
conditions, and they are related to the currents that flow at the crystal surface. These 
conditions define an effective refractive index which is used to study the propagation of 
the ALWs in the material. This index selects which one of the light waves propagates at 
every frequency. 

Here we discuss the effective index obtained from the Pekar's boundary condition. 
This condition is P ex = 0, and it represents the continuity of the excitonic polarization at 
the crystal surface. If there are two light waves the effective refractive index is defined by 
the relation 

ni n 2 

where Q = —Ei/E 2 , and E\ and E 2 are the complex amplitudes of the two waves. From 
Pekar's ABC it follows that Q = (e — n\)/(e Q — n 2 ,) and 

e + n in2 

n cS = ; (10) 

ni + n 2 
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By changing the ABC the behavior of the ratio Q and of the effective refractive index 
changes. Indeed a different choice of the ABC leads to a different choice of E% and E 2 , and 
it affects the propagation. There are many controversial issues related to the consistency 
of the ABCs [T2HT5] . In the rest of the paper we assume that in our case it is possible 
to define n e ff as in ( TTUl) . The effective index n e ff determines which of the two waves is 
propagating or if both of them propagate as follows. 

If the effective index n e ff corresponds to one light wave at lower frequencies and to the 
other at higher frequencies, there is a frequency regime in which n e ff switches from one 
index to the other. This means that both the waves are propagating in that frequency 
range. In this case the KK relations are not satisfied separately on the two solutions 
of ([8]). The standard optics cannot be applied and the correct description is furnished 
by Pekar's theory. The effective index, indeed, satisfies the KK relations and solves the 
apparent violation of causality. 

Standard optics is recovered if the lifetime of the exciton is too short. Indeed in this 
case the exciton does not live enough to modify the linear response of the system and 
a single wave propagates. This possibility is theoretically predicted in Pekar's theory 
because the lifetime of the exciton resonance varies with the temperature. If the temper- 
ature increases up to a critical value T c , also the inverse lifetime increases, and it reaches 
a critical value r c . Below T c the Pekar's theory must be applied, and there is a range of 
frequencies in which both the waves are propagating. Above T c the standard optics holds, 
and the effective index n e ff coincides with only one of the two solutions. 



3 ALW from hydrodynamics 

In this section we observe that the propagation of an ALW is possible in media described 
by relativistic hydrodynamics if the transverse current, that couples to the external EM 
field, has a retarded correlator dominated by a diffusive pole. 

Let us consider a system with a U(l) conserved current J which has a diffusive behavior 
in its transverse part Jt'- {d t — W 2 )Jt = 0, where D is the diffusion coefficient. The 
retarded Green function of Jt has the pole 

with B real, up to higher order in u, k 2 . The analysis of the hydrodynamic equations for 
a relativistic system at finite temperature and chemical potential [T^ITT] shows that the 
transverse current indeed satisfies a diffusion equation; the constants B and V are related 
to the values of the transport coefficients and thermodynamical quantities by 

e is the energy density, p the charge density, P the pressure and rj the shear viscosity. 
If the retarded correlator of the current is dominated by the pole (TTTj) . the transverse 
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permittivity is 



€t(uj, k) = 1 — i 



47rq 2 B 



(13) 



— T>k 2 ) 



where q is the EM coupling. By solving the dispersion equation for (|T3|) we find two 
different solutions for n 2 



There are two different refractive indexes associated to two different light waves. The 
appearance of multiple solutions to the dispersion relation may give origin to the prop- 
agation of ALWs, like in Pekar's theory. The computation of the transport coefficients 
is crucial for the study of the existence and of the propagation of the ALWs in a system 
described by relativistic hydrodynamics. 

4 ALW in holographic optics 

As we discussed above the analysis of ALWs in hydrodynamics requires the knowledge 
of the transport coefficients. There is a class of strongly coupled media in which they 
are exactly calculable, because the Green function follows from the application of the 
holographic principle in string theory. Indeed the holographic correspondence in string 
theory implies that a strongly coupled quantum field theory in <i-dimensional Minkowski 
space-time is equivalent to classical gravity in AdS^+i- In term of this gauge/gravity 
duality the order parameter of the d- dimensional field theory are related to the order 
parameter of a black hole solution in AdS^+i. In this paper we focus on the case d = 4. 
In the hydrodynamical regime we can analitically compute the EM response functions of 
this strongly coupled system, at finite temperature and charge density. If these response 
functions have one or more poles in the (k, u>) plane we can observe the generation of one 
or more ALWs. 

In crystal optics the ALWs are the consequence of the strong spatial dispersion, which 
is associated to the excitonic resonances. The interaction of the excitons with the photon 
is described at macroscopic level by the presence of particular poles in the electromagnetic 
response function e(u,k). In the holographic correspondence we know that the poles of 
the response functions are associated to quasi normal frequencies in AdS gravity [25] . 

Quasi normal modes are solutions to the linearized equations of the classical fluctu- 
ations of the black hole with ingoing boundary conditions. Ingoing boundary conditions 
are chosen because classically the horizon does not emit radiation. The classical fluctu- 
ations are interpreted as small deviations from thermodynamic equilibrium in the dual 
field theory, and the incoming boundary conditions correspond to dissipation. The dis- 
persion relations of this dissipative process are encoded in the poles of the retarded Green 
function, which are located at the (complex) eigenfrequencies of the gravitational quasi 
normal modes. 





(14) 



2V 



7 



We conclude that the quasi normal frequencies in AdS gravity are the holographic 
duals of the excitonic resonances in the strongly coupled field theory. 

In this section we study the propagation of light waves both for transverse and for 
longitudinal waves in media described by a dual gravitational background, in the hydro- 
dynamical regime 0. In the simplest case, at zero charge density but at finite temperature, 
the transverse propagator does not posses any pole, and there is only one propagating 
light wave in the medium. At finite charge density a diffusive pole is generated, and 
indeed this corresponds to a new branch for the refractive index n 2 , which is associated 
to the propagation of an ALW. 

In the longitudinal case, at zero charge density but finite temperature, the correlator 
has a single diffusive pole, and there is a single light wave. Giving a non zero value to the 
chemical potential produces another pole: the second sound. The emergence of this new 
quasi-normal mode generates the second solution of n 2 which is associated to an ALW. 

4.1 Transverse waves without chemical potential: single wave 



Re(n) 

10 r 



Im(n) 

10 r 



CO 



CO 



Figure 1: Re{n) for transverse waves Figure 2: Im{n) for transverse waves 

with T — 1, n — and q — 0.1. with T — 1, /i — and q — 0.1. 

We concentrate on the properties of a medium described by its 3+1 dimensional energy 
momentum tensor T^ v and its 3+1 dimensional U(l) conserved current J M . This subsector 
of the 3+1 dimensional medium is encoded in a 5 dimensional metric g mn and a five 
dimensional gauge field A m , with the action: 

S = ^ J dSF? ( R -^)-^J d 5 x^F mn F™ (15) 

where e is the five dimensional EM coupling constant and 3/1 2 is the cosmological constant. 
We consider the 3+1 dimensional plasma in an homogeneous and isotropic state that is 
in thermal equilibrium, with zero charge density. This ground state of the system is 

4 Observe that the calculation is valid if \k\ — \n\uj is still hydrodynamical 
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described by a particular solution of ( TT5|) . an uncharged AdS black hole: 



ds 2 = -— (dx 2 + dy 2 + dz 2 - f(u)dt 2 ) + -——du 2 (16) 
Ab 2 u v v ' Au 2 f(u) 

where b = an d f(u) = 1 — u 2 . The black hole horizon is at u = 1 and the AdS 
boundary is at u = 0. To compute the linear response of the system at an external 
electric perturbation we need to linearize the equation of motion for the component of 
the gauge field A m (u)e~ wt+lkz parallel to the wave vector k. We need then to solve these 
equations with the infalling boundary condition at u = 1 and arbitrary Dirichlet boundary 
conditions at u = 0. 

The retarded correlator for the transverse current has been computed in [18J at small 
frequencies and wave vector: 

Gr( U ,k) = ~ (17) 
with b = ^f- In this case transverse permittivity €t(u, k) is 

e T {0J, k) = e T (oj) = 1 + 



bco 



The holographic calculation leaves an ambiguity in the definition of the correlator: a 
term proportional to (u 2 — k 2 ). The value of this term has to be fixed by the holographic 
renormalization procedure [19J. The physical requirement we use is that at large oj the 
permittivity goes to its vacuum value: e(co — > oo) = er(uj — > oo,0) — > 1. We computed 
numerically the large u behavior of the correlators and we fixed in this way the value of 
the renormalization constant. Anyway, we observe that all our results are qualitatively 
independent of this renormalization constant. 

Spatial dispersion in ( Tl8|l does not give any effects, and a single wave is expected. 
Indeed the squared refractive index is given by 

n 2 = 1 + (19) 
bus 

and the refractive index, which real and imaginary part are plotted in Figures [1] and [2} is 



V 2b 2 u 2 y/2bcu 

As explained in the Appendix [Bj in this case the sign of the refractive index can be 
studied by comparing Re[n(w)] with the Poynting vector. As claimed in [3] the light wave 
is positively refracted. 

4.2 Transverse waves with chemical potential: ALW 

In this section we study a 3+1 dimensional medium at finite charge density. This ground 
state of the system is described by a different solution of (|T5l) . a charged AdS-RN black 
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Re(n) 



Figure 3: Rein) for transverse waves 
with fi/T = 10 and q = 0.1. 



Im(n) 




Figure 4: Im(n) for transverse waves 
with fi/T = 10 and q = 0.1 



hole, with metric and gauge field given by 

(2 - afl 1 



ds 1 



At — — 



16 b 2 u 



{dx 2 + dy 2 + dz 2 - f{u)dt 2 ) + 



I 2 du 2 
4 u 2 f{u) 



u 
2b 




where 



a + fj, 



f(u) = (1 -u){l + u~au 2 ) 



(21) 



(22) 



As before u = and u — 1 correspond respectively to the AdS boundary and to the black 
hole horizon. The temperature T and the chemical potential \x are given in terms of the 
parameters a and b as 

2 - a [3~ 1 

^=\/n°^ ( 23 ) 



T 
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The retarded correlator for the transverse current has been computed [20] at small fre- 
quencies and wave vector: 



3 a 



ILO 



(2 — a)iu 

4(1 + a)b 2 \-Dk 2 + iu J ~ 8(1 + a) 2 b 



(24) 



2-q 
4ttT 



where D 

This case corresponds to the one studied in [3]. In that case we expanded for small 
k and we ignored strong effects coming from spatial dispersion. In that approximation 
the region of very small frequency cannot be explored, and we found only one wave (with 
negative refraction) propagating in the system. In order to study the effects of strong 
spatial dispersion, we should not expand in the momentum. The transverse permittivity 
is then given by 



e T (co, k) = 1 + 



Ainq 2 



3a 



4(1 + a)b 2 \Dk 2 



iuj 



+ 



(2 



'l + a) 2 b 



(25) 
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whose pole structure is analogous to the one studied in section [3j Indeed, by solving the 
dispersion relation e-r = n 2 , we obtain two different solutions for n 2 . One of the solution 
corresponds to the one studied in [3j. The other instead is the new solution (the ALW) 
associated to the strong spatial dispersion effects. We plot in blue and in red the refractive 
index of the two waves: in Figure |3] the real part while in Figure H] its imaginary part. 

To conclude this section we comment on the sign of the refractive index for the light 
waves. Both light waves are strongly dissipative, as can be seen in Figure HJ The wave 
corresponding to the red line in the Figure [3] is the one studied in [3] . There we studied 
a frequency window where the spatial dispersion effects are small. The Poynting vector 
can then be easily determined and we were able to prove that there is negative refraction. 
The other wave (blu line in Figure |3]) is everywhere in frequencies characterized by strong 
spatial dispersion. Hence the flux energy direction cannot be easily determined and as a 
consequence we need additional informations on the microscopic model to study the sign 
of the refraction (see appendix IB~I). 

4.3 Longitudinal waves without chemical potential: single wave 



Re(n) Im(n) 




CD 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 5: Re(n) for longitudinal waves Figure 6: Im(n) for longitudinal waves 

with T — 1, fj, — and q=0.1. with T — 1, fj, — and q=0.1. 

At finite temperature light waves propagate also in the longitudinal channel. Usu- 
ally longitudinal waves do not propagate in the vacuum because the dispersion relation, 
€l(u) = 0, does not generate any refractive index. Nevertheless at finite temperature the 
background studied above gives origin to a diffusive pole in the current correlator, even 
in absence of finite charge density. This pole in the {k,u) plane implies that a solution 
to the dispersion relation exists and there is a propagating light wave in the longitudinal 
channel. 

The retarded correlator for the longitudinal current has been computed in [TSj at small 
frequencies and wave vector: 

G ^ = ik A b (26) 
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In this case the longitudinal permittivity cl(u;, is 



e L (u,k) = 1 - 27Tq \ (27) 
e z b(tu — k z b) 

The holographic calculation leaves an ambiguity in the definition of the correlator: a 
term proportional to u 2 . The value of this term has to be fixed by the holographic 
renormalization procedure [19]. We use the same requirement of the previous section: 
e(u — > oo) = e L {oj — > oo,0) — > 1. The large u behavior of the correlators has been 
computed numerically, and this fixes the value of the renormalization constant. Once 
again we observe that our results do not depends qualitatively on this renormalization 
constant. The refractive index is computed by solving the dispersion relation ex, (a;, k) = 0. 
There is only a single wave propagating in the medium. In Figure [5] and Figure [6] we 
plotted the real and imaginary part of the refractive index. 

4.4 Longitudinal waves with chemical potential: ALW 



Re(n) Im(n) 




Figure 7: Re(n) and Im{n) for longitudinal waves, with fi/T = 1 (above) and /i/T = 3 
(below) and q = 0.1. 

Here we proceed by analyzing the optics of the system when we add a chemical po- 
tential. At finite charge density there is another pole, corresponding to a new collective 
mode. This mode, called the second sound, is associated to another excitonic resonance. 
Both resonances are peaked around the same frequency u ~ 0, and in this case there are 



12 



two different light waves, two mixing states of the photon with the excitons, propagating 
in the medium. 

The retarded correlator for the longitudinal current has been computed [21 J at small 
frequencies and wave vector: 



G zz (u,k) 
with 



w 2 ( 9a 2(2-a) 2 6 



4(1 + a)b 2 \k 2 - 3lo 2 2(2 + a)bk 2 - 4z(l + a)u - (2 - a) 2 bD m u 



2 

D m = - 



(28) 

(4(1 + a) 3 (l + 4a) 2 log(2-o)-27(2 - a)a 2 (l + 4a) + 2(1 + a) 2 vT+45(2 + o(41o-2)) log (§±^g§)) 



(2-a) 4 (l + 4a) 2 

(29) 

The sound pole (k 2 — 3w 2 ) leads to a second solution for the dispersion relation and hence 
to a second light wave. We fix T = 1 and vary the ratio fi/T, from low to larger values. 
We plot the real and the imaginary part of the refractive index in Figured for \xjT = 1 
and for /j,/T = 3. From the figure we observe that the two solutions n\ and n 2 intersect 
at small frequency for small chemical potential. They do not intersect anymore for larger 
values of //. This effects is associated to the transition between optical regimes, where the 
number of propagating waves changes. We give a more detailed analysis in next section. 

5 Propagation of the ALW in our media 

In this section we discuss the propagation of the ALWs in our media. At finite temperature 
and charge density we have found more solutions to the dispersion relation, leading to 
ALWs. The propagation of ALWs is related to the boundary condition, as explained in 
section [2J 

An important difference between the situation that we studied here and the case of 
crystals, for which the Pekar's' theory has been developed, is that, here, we deal with an 
infinite medium. In crystal optics the ABCs are imposed at the surface of the medium. 
Here we assume that similar conditions can be imposed at infinity such that it is possible 
to define n e fj as in ffTUl) . with e = 1. 

We plot the effective index n e f / for the transverse waves in Figure [SJ From the figure 
we observe that n e ff switches from n\ to n 2 at small frequencies, near the resonance. By 
borrowing the crystal optics physical interpretation, the behavior of n e ff suggests that 
both the light waves propagate. We can imagine an alternative check of this result 
by looking at the propagation of each light wave with respect to the dissipation. This 
property is encoded in the ratio Re(n)/Im(n). We plotted this ratio in Figure ITU1 At 
low frequency both the waves are strongly dissipative, and the ratio Re(n) / Im(n) is 
comparable for them. As the frequency increases the wave selected by n e ff is preferred 
because the absorption is lower for that one. 

An analogous phenomenon is found in the longitudinal waves. We plot in Figure [9] the 
effective index and in Figure [TT] the ratio Re(n) / Im(n) . 
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Figure 8: n e ff for the transverse waves 
at fi/T = 10 and q = 0.1. 
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Figure 10: Ratio Re(n)/Im(n) for the 
transverse waves at fi/T = 10 and q = 0.1. 
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Figure 9: n e ff for the longitudinal waves at 
H/T = 5 and q = 0.1. 
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Figure 11: Ratio Re(n) / Im(n) for the longi- 
tudinal waves at fi/T = 5 and q — 0.1. 



We conclude that both in the transverse and in the longitudinal channel we have 
propagating ALWs. 

5.1 Transitions between different optical regimes 

As explained in section |2j in Pekar's theory for crystals, a transition from propagating 
ALWs to standard optics takes places at some critical value of the order parameter, where 
the lifetime of the exciton become short. 

Here, in the transverse case, at finite temperature and charge density, this transition 
does not occur. At every value of the order parameters there is a range of frequency where 
both waves propagate. Indeed n e ff always switches from ri\ to n 2 for every temperature 
and chemical potential. 

In the longitudinal case the transition to standard optics cannot exist, because, with- 
out taking into account spatial dispersion, there are no light waves propagating in the 
longitudinal channel. Nevertheless we observe here a transition between different optical 
regimes, as can bee seen in Figures WI\ For low values of the chemical potential, n e ff 
coincides with a single light wave in the whole hydrodynamical regime. In this case this 
is the only light wave propagating in the medium. If the chemical potential increases it 
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reaches a critical value \i C rit- For fi > fi crit the effective index n e // switches from n\ to ri2- 
In the frequency range where the transition of n e ff takes place there are two propagating 
light waves. 

The physical explanation of the transition is that the chemical potential is an order 
parameter for the mixing of the two excitons with the photon. For small values of fi 
there is only one light propagating in the medium, and there is no large mixing between 
the two excitons. The single propagating wave is associated to the sound pole in fl28l) . 
which is not dissipative and it is the dominant one at low /i. For larger [i there is larger 
mixing between the two quasi-normal modes, and both can propagate, depending on the 
frequency. 




Figure 12: Optical transition for T = 1 and q = 0.1, at \i = iia-a f or longitudinal waves: 
for ji < n c only one light wave is propagating, while for \i > // c both waves propagate at 
low frequency. 



Conclusions 

In this paper we observed that theories described by relativistic hydrodynamics with a 
diffusive behavior for the currents can give origin to the phenomenon of additional light 
waves. This property has been theoretically predicted by Pekar in crystal optics and then 
experimentally verified. As in the case of crystals, here the ALWs are connected with 
the strong effects of spatial dispersion, i.e. the k dependence of the poles of the Green 
functions. 

We verified our conjecture on the propagation of ALWs in media described by relativis- 
tic hydrodynamics by studying a strongly coupled field theory that can be holographically 
described. In this case the Green functions and the transport coefficients are exactly cal- 
culable thanks to the gauge/gravity duality. 

Note that with our analysis we can only capture the hydrodynamical regime of the 
theory. Indeed at larger frequency and momentum there are other quasi-normal modes. 
They were first computed in [22] for black holes in AdS. These modes exist only for some 
discrete complex frequencies associated to new poles of the Green functions [2"5H2"5] . For 
the transverse and longitudinal R-current correlators the quasi-normal modes have been 
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studied in [26H2E] S In the dual field theory these poles correspond to the presence of 
excitons that can lead to other ALWs, propagating at the various resonance frequencies. 
It would be interesting to perform a detailed study of the existence and the propagation 
of the ALWs out of the hydrodynamical regime. For this computation the use of the 
gauge/gravity correspondence will be fundamental. 

Another issue is the study of the transition to standard optics. In crystal optics it is 
driven by the iT term appearing in the pole of the dielectric tensor. In our case this term 
is set to zero in an ideal system, but it arises if impurities are taken into account. There 
is a general strategy to compute the contribution of the impurities to the Green functions 
in the gauge/gravity duality [30]. Once this contribution is added we expect that there is 
some critical value of the order parameters at which the exciton lifetime is too short for 
the ALWs to propagate. It would be nice to check this prediction in explicit models. 

This paper is a second step in the understanding of optical properties of strongly 
coupled media via holography. Indeed as in [3] the gauge/gravity duality is an useful 
laboratory to verify some general statements common to many systems^. With this mo- 
tivation many other optical properties of strongly coupled systems can be studied. 
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A 6t and ei 

In this Appendix we derive the expression for ej- and in term of the transverse and 
longitudinal correlators. The relevant equation is 

d t E + 4vr J = d t D (30) 

where Di = e^u, k)Ej. 

The vector potential is A4 = P^Ai + PhAi = Af + A\ . The field E t is expressed in 
term of the vector potential A as E^ = F 0i = —iuAi — ikA . From the last equality in 
( 130]) we have 

AnJi = dtDi-dtE, =-iu(D i -E i ) 

= -iu(eij - Sij)Ej = -u(eij - 5ij)(kjA + uAj) (31) 

5 See also for the computation of the residues. 

6 See also [31] . in which the authors studied the optical properties of an holographic superconductor. 
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where hi = (0,0, k). In the linear response theory the Green function is related to the 
current through the relation 



J% = q 2 k)A j + G l0 (u, k)A°) = q 2 (G^co, k)A, - G l0 {u, k)A ) (32) 

We now make use of the identities G z0 = —^G zz and G x q = G y o = G zx = G zy = 0. From 
these relations and from ( 13~T|) and (|32|) we obtain 

AlT 

€ij(uj, k) = 8 i:j - — q 2 Gij(u, k) (33) 

By projecting this relation on the transverse and longitudinal channels we obtain 

Air Air 
e T (u,k) = l- — q 2 G T (u,k) , e L (u,k) = l- — q 2 G L (u,k) (34) 



B Poynting vector 

In the paper we showed different physical systems exhibiting the phenomenon of ALW. It 
is interesting to determine if a propagating wave shows negative refraction too. There is 
negative refraction when the phase vector is opposed to the energy flux vector. The sign 
of the phase velocity is given by Re[n(u;)] while the energy flux is given by the Poynting 
vector. However, the Poynting vector is not always well defined [Sj. 

In the case of small dissipation, even if spatial dispersion is present, the Poynting 
vector is oriented as the group velocity v g , and has the simple expression 

5 = Re( E *AB-|^E*E,) (35) 

Negative refraction in this context has been studied in [7]. 

When dissipation effects are not negligible, the Poynting vector is not in general di- 
rected as the group velocity. If spatial dispersion effects are small, we can compute it 
using the classical e, \x approach \j (see [32|[33]). 

S = Re \E\ 2 (36) 

In [3] we showed that this expression is equivalent to (|35|) for transverse waves, by ex- 
panding the permittivity at second order in the wave vector. In that case we can studied 
negative refraction by comparing the sign of the Poynting vector with the sign of Re[n(o;)]. 

In the general case of strong spatial dispersion and strong dissipation we cannot iden- 
tify precisely the direction of the flux of energy [8]. Hence we cannot determine if there 
is negative refraction. Except for the wave also analyzed in [3], all the other waves we 
studied in this paper fall in this category. 

7 Where [i here refers to the magnetic permeability. 
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